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Introduction 

There is a strong similarity between the expansions of the heat kernel as worked 
. out by people in Riemannian geometry in the seventies (starting with the famous 

^ ! "Can one hear the shape of a drum" by Mark Kac |Kacj . the Berger paper Berger 



> 



and the Mc-Kean-Singer paper |McK-Si] ) and the so-called semi-classical trace 
formulas developed by people in semi-classical analysis (starting with Helffer- 
Robert |He-Roj ) . In fact, this is not only a similarity, but, as we will prove, 
each of these expansions, even if they differ when expressed numerically for some 
^ ! example, can be deduced from the other one as formal expressions of the fields. 

t:;^- I Let us look first at the heat expansion on a smooth closed Riemannian mani- 

fold of dimension d, {X,g), with the (negative) Laplacian aJ^. The heat kernel 
O ' e{t,x,y), with t > and x,y G X, is the Schwartz kernel of exp(tAg): the 

solution of the heat equation Ut — AgU = with initial datum uq is given by 

; u{t,x)= e{t,x,y)uo{y)\dy\g . 

The function e{t,x,x) admits, as t — )■ 0"'', the following asymptotic expansion: 

e(t, X, x) ~ {Anty^/^ (l + ai{x)t + ■■■ + ai{x)t^ + ■ ■ ■ ) . 

The a/'s are given explicitly in |Gilkey2] , page 201, for / < 3 and are known for 
/ < 5 [Avramidil IvdV] . See also the related works by Hitrik and Polterovich 
[Hit IHi-Po-ll IHi-Po-2l IPoj . They are universal polynomials in the components of 
the curvature tensor and its co- variant derivatives. For example oq = 1, oi = Tg/6 
where Tg is the scalar curvature. 



*Institut Fourier, Unite mixte de recherche CNRS-UJF 5582, BP 74, 38402-Saint Martin 
d'Heres Cedex (France); yves.cohn-de-verdiere@ujf-grenoble.fr 

""^In this note, we wiU not foUow the usual sign convention of geometers, but the convention 
of analysts 
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The previous asymptotic expansion gives the expansion of the trace by inte- 
gration over X and has been used as an important tool in spectral geometry: 



trace(e*'^'') = / e{t, x, x)\dx\g = >^ 



CO 



where — Ai = < — A2 < ■ ■ ■ < —Xk < ■ ■ ■ is the sequence of eigenvalues of — A^, 
with the usual convention about multiplicities. If c? = 2, this gives 

trace(e*^'') = ( Area(X) + + 0{t^)] , 

Ant \ 6 / 

where x{^) is the Euler characteristic of X. 

There is an extension of the previous expansion in the case of Laplace type 
operators on fiber bundles: the coefficients of the expansion are then polynomials 
in the co- variant derivatives of the curvature of the metric and of the connection 
on the fiber bundle. The heat expansion can be re-interpreted as an expansion 
of the Schwartz kernel of f{—h^Ag) on the diagonal x = y in powers of h with 
f{u) = exp{—u) and t = h"^. This is a particular case of the semi-classical trace. 

Let us describe the semi-classical setting in the fiat case: H;i is a self-adjoint 
/i— pseudo-differential operator with Weyl symbol H{x,C,) in some open domain 
X in M'^, or more generally on a Riemannian manifold. Let / G iS(M) and look 
at f{Hti). Under some suitable assumptions (ellipticity at infinity in ^) on H, 
f{Hfi) is a pseudo-differential operator whose Weyl symbol f*{H) is a formal 
power series in h, given, using the Moyal product denoted by ^, by the following 
formula (see jGracia] for explicit formulas and Section 4.2 for a proof, see also 
[Charles] ) at the point Zq G T*X: 



rmzo) = {2nh)-'^ [y, \/\h{z,)) {h - H{z,)r ' (^o) j . (i) 

From the previous formula, we see that the symbol of f{Hn) at the point z 
depends only of the Taylor expansions of H at the point z and of / at the point 
H[z). In the paper |H-P] . the authors study the case of the magnetic Schrodinger 
operator whose Weyl symbol is Hay{x,^) = X]j=i(0 ^ %(^))^ + and show 

that the Schwartz kernel of f{Hfi^ay) at the point (x, x) admits an asymptotic 
expansion of the form 



[f{Hn,ay)]{x,x) = {2nh)- 



00 I r 



-d 

^ j=Q \ 1=0 

where the Q'^'^ {x, .^)'s are polynomials in calculated from the Taylor expansions 
of the magnetic field B = da and V at the point x. The proof in |H-Pj uses a 
pseudo-differential calculus adapted to the magnetic field. 
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We will give a simplified version of the expansion replacing the (non unique) 
Qj'/^(x,^)'s by functions Pj^'^(x)'s which are uniquely defined and are given by 
universal 0((i)-invariant polynomials of the Taylor expansions of B and V at the 
point X. We present then two ways to compute the Pj { 's: 

• we can first use Weyl's invariant theory (see Gilkey's book |Gilkey2| ) in 
order to reduce the problem to the determination of a finite number of 
numerical coefficients; then simple examples, like harmonic oscillator and 
constant magnetic field, allow to determine (part of) these coefficients. 

• The P^i^'s are related in a very simple way to the coefficients of the heat 
expansion; it is possible to compute the P- { 's from the knowledge of the 
a^'s for j + 1 < / < 3j. This is enough to re-compute the coefficient of 
and also, in principle, the coefficients of h'^ in the expansion, because the 
a^'s are known up to / = 6 in the case of a fiat metric (see |vd V] ) . 

In this note, we will first describe precisely the semi-classical expansion for 
Schrodinger operators (in the case of an Euclidian metric) and the properties 
of the functions Pj^'^(x)'s. Then, we will show how to compute the Pj^'^(x)'s 
using an adaptation of the method used for the heat kernel (Weyl's Theorem on 
invariants and explicit examples). Finally, we will explain how the a^'s are related 
to the Pj^^^ixys. This gives us two proofs of the main formula given by Helffer 
and Purice in |H-Pj : this paper was the initial motivation to this work. 
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the statements and proofs more precise. Many thanks also to Johannes Sjdstrand 
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me to present them in an Appendix to this paper. 

1 Semi-classical trace for Schrodinger operators 

In what follows, X is an open domain in M'^, equipped with the canonical Eu- 
clidian metric, and Q^{X) will denote the space of smooth exterior differential 
forms in X. Let us give a Schrodinger operator, with a smooth magnetic field 
B = ^i<j<j<d bijdxi A dxj (a closed real 2-form) and a smooth electric potential 
V (a real valued smooth function) in X. We assume that V is bounded from be- 
low. We will assume also that the 2-form B is exact and can be written B = da 
and we introduce the Schrodinger operator defined by 

i^M,v = $^(-^-%(x)j +V{x) . 
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The Weyl symbol of Hn,ay is Hay{x,i) = ||^ - a{x)f + V{x). We denote by 
Hri,a,v a self-adjoint extension of Hr^ay in L'^i.^i \dx\). Let us give / G S{M) and 
G C^(X) and consider the trace of (j)f{H^^ay) as a distribution on X x M (the 
density of states): 

Trace (0/(ifr*,a,y)) = / ^r»,ay(fi')(a;)0(x)|cix| , 

where Zn^ay{,g){x) is the value at the point (x, x) of the Schwartz kernel of 

Theorem 1 We have the following asymptotic expansion in power of h: 

l=3j 



Zn,ay{g){x) ~ {2nhy 
We have the explicit formulas 



[ /(lief + n^))Mei + E^'M E 



By, 

j+i 



/ f^\ur+v{x)m\ 



Pi,^ = -\{^v+\\Br) 

P^/ = -Y^\\^v\\\ 



= (8||V5f + \\d*Bf + 12{AB\B) + 3AV) . 

Here ||-B|p = J2i<i<j<d^ij' ^* '■ ^^(-^) fi^(X) the formal adjoint of d used 
in the definition of the Hodge Laplacian on exterior forms. Ifd = 3, \\B\\ is the 
Euclidean norm of the vector field associated to B. 

The P^i^ (x) are polynomials of the derivatives of B and V at the point x. 
Moreover, if X, fi, c are constants and we define X*{f){x) = f{Xx), we have the 
following scaling properties: 

1. P^f^^^'^*^^\x) = X^^Pff{\x). This will he used with x = 0. 

By+C, N _ ryBy, 



3- P-r^'ix)=P-{ (x 

^B 
hi 



I p-ri^=p'ri^- 



5. ThePlff 's are invariant by the natural action of the orthogonal group 0{d) 
on the Taylor expansions of B and V at the point x. 

Remark 1 From the statement of the previous Theorem, we see that the expan- 
sion of the density of states is independent of the chosen self-adjoint extension. 
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As a consequence, we can get the following full trace expansion under some 
more assumptions: 

Corollary 1 Let us assume that Eq = iniV < Eoo = liminix^ax V{x) and that 
we have chosen the Dirichlet boundary conditions. Let f G C^{] — oo, Eoo[) , then 
the trace of f{H^ a,v) admits the following asymptotic expansion 

Traceifilh^)) ^ {27rhr' [ ([ f + V {x)M\ + ■ ■ ■ 

oo l=3j 
i=l l=j+l 

The coefficient of fi^ can he written as 

-i- / f^'\Ur + V{x)) {AV{x) + 2\\B{x)r) \dxd^\ . 
The expansion follows from [He-Ro] . An integration by part in x gives 

[ f^'\ur+v{xm'^v{x)r\dx\ = - [ f^'\ur+v{x))Av{x)\dx 

J X J X 

2 Existence of the /i— expansion of Zfi^ay 



Using Theorem [2] in the Appendix, we can work in R'' with a and V compactly 
supported. The existence of the expansion is known in general from |He-Ro] and 
the calculus of the symbol of f{Hfi^a,v)- We get 




where the Qjj{x, ^)'s are polynomials in the Taylor expansion of Hay at the point 
(x, ^). The previous expansion is valid for any (admissible) pseudo-differential 
operator. In the case of Schrodinger operators we can make integrations by part 
in the integrals / f^''\Hay{x,^))Qj^i{x,^)\d$,\ which reduces to a similar formula 
where we can replace the Qjj{x,^)^s by Pj^;(x)'s. This is based on the expansion 
of Qjj as a polynomial in ^ in powers of — a): odd powers give and even 
powers can be reduced using 

4 ((e, - a,)f^'\Hay)Lid^^)d^) = 2||e, - a,ff^'+'\Hay)d^ + f^'\Hay)di . 

We have only to check that the powers of ^ in Qj^i{x,^) are less than /: this 
is based on Equation ([1]). The coefficients of the /— th Moyal power of Hay{z) — 
Hay{zQ) are homogeneous polynomials of degree I in the derivatives of Hay{z). 
At the point z = Zq only derivatives of order > 1 are involved. They are all 
of degree < 1 in Using Gauge invariance at the point x (Section [3]), we can 
assume that a{x) = 0. 



f\ur+v{x)m\ I \dx\ 
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3 Gauge invar iance 



If S* : X — J- M is a smooth function, we have 
and 

Hence, we can chose any local gauge a in order to compute the expansion: 
using the synchronous gauge (see Section Hj), we get the individual terms 

for the expansion, where the P^i^ (x) depend only of the Taylor expansions of B 
and V at the point x. 



4 The synchronous gauge 

The main idea is to find an appropriate gauge a adapted to the point Xq where 
we want to make the symbolic computation. In a geometric language, we use 
the trivialization of the bundle by parallel transportation along the rays: the 
potential a vanishes on the radial vector field H. Here, this is simply the fact that, 
for any closed 2-form B on M?, there exists an unique 1-form a = Yl'j=i djdxj so 

that da = B and ^^=1 Xjaj = 0. 

We will do that for the Taylor expansions degree by degree. In what follows 
we will use a decomposition for 1-forms, but it works also for /c— forms. 

Let us denote by the finite dimensional vector space of /c-differential forms 
on whose coefficients are homogeneous polynomials of degree N and by = 
Yl'j=i ^iaf" the radial vector field. The exterior differential induces a linear map 

from fi^ into and the inner product i{W) a map from into They 

define complexes which are exact except a.t k = N = 0. Moreover, we have a 
situation similar to Hodge theory: 

This is due to Cartan's formula: the Lie derivative of a form u G Q% satisfies, 
from the direct calculation, Cyyu = {k + N)u, and, by Cartan's formula, = 
d{L{W)u) + L{W)duj. So 



^This gauge is sometimes called the Fock-Schwinger gauge; in |A-B-P| . it is called the 
synchronous framing 
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It remains to show that this is a direct sum: ii u = da = L{W)'y, we have 
l{W)u = and du = 0; from the previous decomposition, we see that u = 0. 
Let us denote by J^u, where a; is a differential form of degree k, the form in fi^ 
which appears in the Taylor expansion of u. 
We get 

Proposition 1 If P{J^a, J^a, ■ ■ ■ , J^a) is a polynomial in the Taylor expansion 
of the 1-form a at some order N which is invariant by a ^ a + dS, P is inde- 
pendent of J^a and 

P {J'a, ■ ■ ■ , J^a) = P (^ljh{W)B, ^±^J^,(W)B^ 

is a polynomial of the Taylor expansion of B to the order N — 1. 

5 Properties of the Pj/s 

5.1 Range of / for j fixed 

From the scaling properties, we deduce that, in a monomial 

D'^'B,,,,, ■ ■ ■ D'"'B,^,,^D^'V ■ ■ ■ D^-V , 

belonging to Pj^, we have k + 2m = 2{l— j) and + |ai| H — ■ + |afc| + + ■ — \- 
|/3m| = 2j. Moreover, for j > 1, /c + m > 1 and > 1. Hence j + 1 < / < 3j. 
The previous bounds are sharp: take the monomials A^V and || VV^H^-' which give 
I = i + 1 and / = 3j. 

5.2 Invar iance properties 

1. Let us assume that we look at the point x = and consider the operator 

= /(a^^;). We have 

The same relation is true for any function f{Hh,A,v) and then we have, 
looking at the Schwartz kernels and using the Jacobian /x'^ of D^: 

p,7io) [ f^^\ur+vmm = fi-'^p^f''^''''^{o) [ f^^\ur+vmm. 

2. We have 
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3. Changing V into V + c gives a translation by c in the function / but does 
not change the P^^'^'s. 

4. Changing B into —5 gives a complex conjugation in the computations. 
The final result is real valued. 

5. Orthogonal invariance is clear: an orthogonal change of coordinates around 
the point x preserves the density of states. 



5.3 The case d = 2 

We deduce from the scaling properties and invariance by the orthogonal group, 
that there exists constants a^, bd, q so that P]^2^(a;) = adAV+bd\\B\\^, Pi 3(0;) = 



6 Explicit examples 

The calculation for the harmonic oscillators and the constant magnetic fields 
allows to determine the constants a^, bd, Cd- 



6.1 Harmonic oscillators 

Let us consider f2 = —h^-^ + x"^ with d = 1. The kernel of P(t, x, y) of exp(— 
is given by the Mehler formula: 

P{t,x,y) = {27rhsmh{2th))-hxp (^--^—L-^ {cosh{2th){x^ + y^) - 2xy)^ . 
Hence 

P(t,x,x) ~ {27ih)''e~''^' (^j e-'^'di^ {l - h^{f - t^x^) /?> + Q{h^)) . 

Hence Pi,2(x) = -V"{x)/Q, Pi,3(x) = -V\xf/12,. 

Similarly, in dimension c? > 1, we get Pi, 2(2^) = — Ay(x)/6, Pi.3(x) = 
-||Vy||Vl2,. 



6.2 Constant magnetic field 

Let us consider the case of a constant magnetic field B in the plane and denote 
by Q{t,x,y) the kernel of exp(— tif^^o)- We have (see [A-H-Sj ) 
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Hence the asymptotic expansion 

Q{t,x,x) = {2TTh)-' J exp{-tUf)\d^\{l-t^h^By6 + 0{h'')) , 

hence Pi,2(a;) = Pi^six) = 0. 

Using the normal form B = bi2dxi A dx2 + b^idxs A dx^ + ■ ■ ■ , we get in 
dimension d > 2, Pi,2{,x) = — ||i?p/6, Pi,3(x) = 0. 



7 Heat expansion from the semi-classical expan- 
sions 

We have tHi a,v = H^^^^ ^y. Using the expansion of Theorem [1] with f{E) = 
, we get easily the point-wise expansion of the heat kernel on the diagonal as 
t 0+: 



[eM-m,a,v)]{x,x)r.^—^^e-^yi^) 



(=0 \V3<j</-l 



In particular, ai{x) = —V{x) and the coefficient a2{x) is given by 
a2{x) = lv{xf-lAV{x)-l\\B{x)r. 

This formula agrees with Equation (3) of Theorem 3.3.1 in |Gilkey2| . 
This gives another way to compute the Pj/s: if, as power series in t, 

oo / oo 

1=0 \l=0 

we have 

J2 Pff{x)=h{x). 
l/'i<j<l-^ 

Pfi^ is the sum of monomials homogeneous of degree 2{l — j) in bi where B and 
its derivatives have weights 1 while V and its derivatives have weights 2. 

The heat coefficients ai on flat spaces are known for / < 6 from |vdV] . This 
is enough to check the term in h"^ (uses 02 and as) in |H-P] and to compute the 
term in h!^ in the semi-classical expansion (uses the a^'s for 3 < / < 6). 

We have also a mixed expansion writing tHn^ay = H^j^ ^y, we get a power 
series expansion in powers of h and t valid in the domain h^t — ?■ and < t < to 
for the point-wise trace of exp{—tH^^a,v)- 

^ ' \ i>i, j+i<«<3j / 
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This shows that the integrals J^V{x)'^\dx\ and P^i^ {x)\dx\ are recoverable 
from the semi-classical spectrum. 

Appendix: functional calculus in domains and 
self-adjoint extensions (after Johannes Sjostrand) 

The content of this Appendix is due to Johannes Sjostrand. I thank him very 
much for this contribution. 

Let X C M'^ be an open set, we say that a linear operator A is a "^DO in X, 
with Weyl symbol a, if, for any compact K C X, A acts on functions supported 
in K ^DO of Weyl symbol a. 

Theorem 2 Let Hu ay be a Schrddinger operator with magnetic field given by 
Hh,ay = Y^{-^-aj{x)\ +V{x) , 

defined in some open domain X d W^. We assume that a and V are smooth 
in X and that V is bounded from below, so that Hu ay admits some self-adjoint 
extensions on the Hilbert space L'^{X, \dx\). One of them will be denoted by Hfi ay. 
Then, for any f G iS(M), f{Hfi^ay), given by the functional calculus, is a semi- 
classical ^DO in X whose symbol is given by Equation (|7]j and is independent 
of the chosen extension. 

The proof uses a multi-commutator method already used by Helffer and 
Sjostrand. Proof. - 

We introduce, for s G M, the semi-classical (^-dependent) Sobolev 
spaces 

m := {u G S'iR") I ||Op,(l + UWYMl^ < oo } 
with the norm 

\\u\U:= \\Op,{l + ufr^ML^ . 

The (^-dependent) norm ||A||si,s2 is the norm of A as linear operator 
from to A linear operator K is smoothing if, for all Si, S2, 
Il-^l|si,s2 ~ 0{h°°) . This implies that the Schwartz kernel of K is 
smooth with all derivatives locally 0{h°°). We have the 

Lemma 1 Let Y be an open set m M.'^. Let Pj = Pj{h), j = 0, 1 
be two self-adjoint operators on Hilbert spaces Hj = L'^{Xj, \dx\) with 
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y CC Xo C Xi C M*^ and with domains Vj so that C^{Y) C Vj C 
Hj. Let us assume that, on C^{Y), Pq = Pi = Hfi^ay{= P)- 

Then, for any f G C^(R), f{Po) — /(-Pi) is smoothing on Y . In 
particular, the densities of states [f (Pj)]{x , x) , j = 0,1, coincide in 
Y modulo 0{h°°). 

Assuming Lemmadl Theorem [2] follows by extending a and V smoothly 
outside Y so that they have compact support in Mf^. We take Y GG 
X = Xf) G M.'^ = Xi. It follows that Pi is essentially self-adjoint and 
the functional calculus for Pi follows then easily from |He-Ro] . The 
result is valid even for / G S(R) because is dense in S and the 
result of |He-Ro] is valid for f G S and the resulting formulas for the 
symbols are continuous w.r. to the topology of S. 

□ 

Proof of LemmaUl- If X ^ C^{Y), then, for z ^ M and j, /c G {0, 1}, we have 
on L^iY): 

{P, - z)-' ox = Xo{Pk- z)-' - {P, - z)-'[P, xm - z)-' (2) 

Let xo < Xi < ■■■ < Xtv with, for / = 0, ■ ■ ■ , iV, x« e C^{Y) and, for / = 
0, ■ ■ ■ , iV - 1, - Xz+i) = 0. By iterating Equation and using Xi+i[P, Xi] = 
[P,Xi], we find: 

(Pi - z)-' oxo = xi° {Po - zy'xo - X2 o (Po - zy^[p, xi]{Po - zy'xo + ■■■ 
±xn{Po - zy'[p,XN^i]{Po - zy'[p,XN^2] ■■■{Po- zy'xo 

T{Pi - zy^[P, xn]{Po - zy ■■■{Po- zy^xo 

Let us give now xo, ^ G C^{Y) with disjoints supports. By choosing the x'l^ 
for / > with supports disjoint from the support of ^/j, we see, using Equation 
(E]), that, for any N, 

||^(Pi-z)"ixo||o,2 = 0(n^|53z|-(^+i)) . 

The standard a priori elliptic estimates 

||'«L+2,ni < C (||(P - z)u\\s,n2 + ll'^IUnJ 

for z E K CC C and Qi CC Q2 CC M'^, allow to prove that, for any N, s, there 
exists M{N, s) so that 

||^(Pl - zy\o\\s,s+N+2 = 0(n^|3z|-^-^(^'^)) (3) 

Let X G C^{Y) so that x = 1 on the support of xo- Let us apply a multiplication 
by Xo to the right and to the left to Equation ([2]) and choose ip with support 
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disjoint of xo so that [-P, x]"^ = Inserting ip this way in Equation ([2]), we 

get, using Equation ([3]): 



and, for any TV, there exists M{N) so that ||/r||_Ar,7v = 0{h^'^z-^^'^^^). We now 
apply the formula (known to some people as the "Helffer-Sjostrand formula", 
proved for example in the book |Di-Sj| , p. 94-95), valid for / G C^(M) and / an 
almost holomorphic extension of /: 



where dL{z) is the canonical Lebesgue measure in the complex plane. From this, 
we see that /(-Po) — /(-Pi) is smoothing in Y. 
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